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DIFFUSION PROCESSES 



ANDREA POSILICANO AND STEFANIA UGOLINI 

Abstract. We suggest a rigorous definition of tlie patliwise flux 
across the boundary of a bounded open set for transient finite en- 
ergy diffusion processes. Tfie expectation of sucfi a flux has the 
property of depending only on the current velocity v, the nonsym- 
metric (with respect to time reversibility) part of the drift. In the 
case where the diffusion has a limiting velocity we define the as- 
ymptotic fiux across subsets of the sphere of radius R, when R 
tends to infinity, and compute its expectation in terms of v. 

Keywords: diffusion processes, (random) flux across surfaces, cur- 
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1. Introduction 

In a previous paper [THj the authors gave a pathwise probabihstic 
versions of the Scattering-into-Cones and Flux-across-Surfaces theo- 
rems in Quantum Mechanics and then recovered the known analytical 
results by taking suitable expectations. 

Here we extend the main probabilistic results contained in to 
a large class of Markovian diffusions with no a priori connection with 
Quantum Mechanics. 

We consider diffusions on M'^ which are weak solutions of s.d.e.'s of 
the form 

dXt = b{t,Xt)dt + dWt, (1.1) 

where Wt is a standard Wiener process and the drift vector field bt{x) = 
b{t, x) satisfies 

VT>0, e[ \\b{t,Xt)fdt<+oo. (1.2) 



^0 

The condition (1.2) is a finite energy condition in the sense of Follmer 
0. By Girsanov theory one proves that (1-2) is equivalent to a finite 
entropy condition: the probabilitity measure describing the weak so- 
lution of (1.1) has finite relative entropy with respect to the Wiener 
measure. By definition of relative entropy, this fact implies absolute 
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continuity and therefore the distribution of X^, t > is absolutely con- 
tinuous (with respect to Lebesgue measure) with some density function 
Pt- 

A very peculiar consequence of (1.2) is that not only the Markovian 
property is preserved under time reversal, but also the diffusion one. 
We will call this property time reversibility of the diffusive character. 
The class of diffusion processes individuated by this invariance prop- 
erty was firstly proposed by Nelson in 1966 ^01 in the framework of 
Stochastic Mechanics. Such diffusions were then studied by Zheng and 
Meyer who called them semimartingales dans les deux directions du 
temps Within this class, Carlen in 1984 |^ (see ^ for an alter- 
native proof) solved the existence problem of weak solutions of (1.1) 
in the case of (unbounded) drift fields satisfying (1.2). Successively, 
FoUmer [H] gave a very elegant characterization based on the relative 
entropy approach. Under time reversal the process solution of (1.1) is 
again a solution of a s.d.e. of type (1.1) with some dual drift field bl 
(and, of course, another standard Wiener process) which satisfies the 
relation: 

bt{x)-bl{x)=V log pt{x) (1.3) 

The duality relation (1.3) allows to introduce a relevant decomposition 
of the drift field as the sum of two vector fields : 

bt + bl bt - bl 

^* = ^-' "* = ^- 

called respectively current and osmotic velocity JT]. In the symmetric 
case, Wf = and thus bt = Ut is of gradient type according to (1.3). 
Therefore the current velocity Vt represents the non symmetric part of 
the drift field. We will see that only the current velocity is involved in 
the expression of the flux across surfaces by the diffusion paths. 

A useful consequence of the time reversibility of the diffusion prop- 
erty is then the validity (in the weak sense) of the continuity equation 
for the couple (p, u): 

^^Pt = -V -{ptVi). (1.4) 

In this paper we give a rigorous definition of the pathwise flux across 
the boundary of a bounded set by a transient Markov diffusion process 
solution of (1.1) and satisfying the energy condition (1.2). 
Given a bounded open set D let us consider the function 

NaD:=N^D-Naj,^ (1.5) 

where Ng^{'y) (resp. Ngj^{'y)) denotes the number of inward (resp. 
outward) crossing by the path 1 1-^ ^4(7) of the boundary dD, 7 being 
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the point in the probabihty space. It is not a local time because we 
need to distinguish outward from inward crossing in order to have the 
pathwise analogoue of a net flux. The problem is that almost surely 
the diffusion Xt intersects dD on a set of times that has no isolated 
points and is uncountable. Therefore the definition of Nqd given above 
makes no sense in general. However, by a suitable redefinition of Nqd 
as the total mass of the almost surely compactly supported random 
distribution —-^XDi^t) , where xd is the characteristic function of 
the set D (see section 3 for details), we can give a rigorous definition of 
the pathwise flux across dD. Then, by using the continuity equation 
(1.4), in the case where dD is a regular hypersurface we can compute 
(see theorem 3.2) the expectation of the pathwise flux Nqd in 
terms of the current velocity v, obtaining 



where n denotes the outward unit normal vector along dD and a is the 
surface measure. We interpret ^qd as the flux of Xt across dD. Note 
that, by our choice of signs in (1.5), D is a source if ^qd > and is a 
sink if < 0. 

Of course the definition of flux given above does not extend to the 
case of a hypersurface which is not a boundary. This restriction can 
be avoided, at least asymptotically, in the case where the diffusion Xt 
has a limiting velocity, i.e. 



exists almost surely for some non zero random vector Voo- A simple but 
general condition giving (1.7), which is again a finite entropy condition, 
was obtained by Carlen [3] (see Theorem 4.3 below). 

Suppose that S is an open subset of the unit sphere with (9S a finite 
union of manifolds. In order to define the asymptotic flux across S 
we consider the cone Cs := {Ax G M'^ : x G S, A > O}. Using (1.7) 
(see section 5 for the details) we can then define N^, the asymptotic 
pathwise flux across S, by the limit t cxd of the mass N^,^ of the 
almost surely compactly supported random distribution ■^xcnB'^iXt) , 
Br being the closed ball of radius R. Note the change of sign in the 
definition of A^s^ with respect to Nqq. This is consistent with the 
fact that the exterior normal to := Cs HS'r, Sr being the sphere of 
radius i?, coincides with the interior normal to the boundary of Cfl-B^. 
No confusion can arise between the two different definitions since 
is never the boundary of an open subset of M''. 




(1.6) 



lim - 

tToo t 



t = V, 



oo 



(1.7) 
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We show (see theorem 5.2) that is well defined since almost surely 
one has 

:= \im N,:, ^ XcM ■ (1-8) 

R]oo 

Moveover, if as before we define the flux by taking the expectation of 
the corresponding pathwise object, we have (see theorem 5.3) 

r+oo n 

$|=lim / / p{t,x)v{t,x) ■n{x)d(T{x)dt . (1.9) 

The conditions required in order to obtain the stated results are, besides 
(1.1) and (1.2), that pt and Vt belong to the Sobolev space H^{M^) in 
order to obtain (1.6), (1.7) to get (1.8), and moreover we assume 

E ! \\Vv{t,Xt)fdt < +00 

for some T > 0, to get (1.9). 

Finally let us remark that our definition of a pathwise flux across 
dD is euristically equivalent to the ill-deflned Stratonovich stochastic 
integral 

(^j 5{Xt-x)n{x)da{x)^odXt. 
Indeed, proceeding euristically, one obtains 

--XD(^t) ^ • ^XD{Xt) dt 

= ^ . n(Xt) 59D(Xt) dt 
^lo (/ ~ ^'^ ^^^-^ da(x)^ o dXt . 



2. The class of finite energy diffusion processes 

Consider the measurable space (f2, ^), with Q — C(M+; M"^), JF the 
Borel (T-algcbra. and let {Q, JF, JF^, Xf) be the evaluation stochastic 
process Xt{'y) :— j{t), j & fl, with J^t — cri^s, < s < t) the natural 
flltration. 

Let us suppose that: 

HI) there exists a Borel probability measure P on (fi, JF) such that: 



- (Q, J^, J^ti ^t, P) is a Markov process; 
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- Wt := Xt — Xq — / b{s, Xs) ds is a (P, JFt)-Wiener process, i.e. P is 
Jo 

a weak solution of the stochastic differential equation 

dXt = b{t,Xt)dt + dWt (2.1) 

with initial distribution /x = P o Xq^, /i <^ A, A the Lebesgue measure 
on R'^; 

H2) the adapted process bt satisfies 

VT>0, e[ \\b{t,Xt)fdt <+oo (2.2) 
Jo 

where E denotes the expectation with respect to the measure P. 
Remark 2.1. Since 

d¥ 



i/^^(P,P^) :=Ellog 



dF^ 



iE^^||6(t,X,)|pdt 



where P^ := J /i(x)Pjf , with denoting the Wiener measure starting 
from X, see jH], the finite energy condition (2.2) is equivalent to a finite 
relative entropy condition. Thus (2.2) implies that P is absolutely 
continuous with respect to P^, and so Xt admits a density function pt 
for any t > 0. 

As a consequence of HI and H2, the Markovian diffusion Xt preserves 
the diffusion property under time reversal. Indeed by FoUmer P one 
has: 

Lemma 2.2. Under the hypothesis HI, H2, defining P := Poi?, where 
R is the pathwise time reversal on C{[0,T]]W^) , R{'~f){t) := 7(T — t), 
there exists an adapted process bt such that: 

Wt:=Xt-Xo- [ b{s,X,)ds (2.3) 
Jo 

is a (F, J-'t) -Wiener process. 

Remark 2.3. Lemma 2.2 states that the finite energy (entropy) con- 
dition (2.2) is a sufficient condition for the time reversibility of the 
diffusion property. The proof is based on the fact that the finite en- 
tropy condition is invariant under time reversal. The extension to the 
infinite dimensional case is in j7j. Sufficient conditions are also given 
in [H]. Sufficient and necessary conditions for reversibility of diffusion 
property, in the case of Lipschitz drift fields, are investigated in jU]. 
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It is well known that the drift field can be seen as a stochastic forward 
derivative in the sense of Nelson ^U], ^1]. In particular from (2.1) and 
(2.2) it follows that (see 0): 

bt = lim ]-E[Xt+h - XtlJ't] in L^F) . 

Analogously, from (2.3) one also has: 

bt = lim ]-t[Xt+h - Xt\Tt] in ^^(P) , 
hio n 

where E denotes the expectation with respect to the measure P. For our 
approach it is convenient to work with the same probability measure 
P as proposed by Nelson ^U]. To this end we write: 

bt = lim —EiXr-t - XT-t-h\:FT-t] o R in ^^(P) , 
hio h 

where J-'t = cr{Xs, s > t) is the natural future filtration. 

Since the Markov property is preserved under time reversal also the 
dual drift is given by some measurable function 64(7) = b{t, {Xt{'j)). 
Let us define 6j (x) = —br-tix) so that, as already obtained in [TT] : 

b; = lim - ]-E[Xt - Xt-h\:Ft] in (P) , 

and the following relation holds: 

bt{x)-b;{x)=V log pt{x) (2.4) 

between the drift field and its dual (see and for the non Markovian 
case jni). 

The duality relation (2.4) allows to introduce the decomposition 

bt = ut + Vt, bl = -ut + Vt 

where 

bt + b1 bt - 6* 

Vt = Ut = ^- 

are called current and osmotic velocity respectively |lUj . 

In the symmetric case, Vt = 0,bt = Ut,b^ = —Ut, thus the drift field 
coincides up to the sign with its dual and, according to (2.4), it is 
of gradient type. Therefore the current velocity Vt represents the not 
symmetric (with respect to time reversal) part of the drift field. We 
will see that only the current velocity is involved in the expression of 
the flux across surfaces. 

An important consequence of the time reversibility of the diffusion 
property is the validity of the continuity equation (in the weak sense) 
for the couple (p, f). 
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Indeed, recalling the Fokker-Planck equation associated with (2.1): 
d 

— pt^-V-{ptbt) + Apt 
and the Fokker-Planck equation associated with (2.3): 

^Pt = -V-(pA*)-Apt 
and putting together the two equations one has: 

^p,^-V-{ptVt) (2.5) 
where the definition of the current velocity has been used. 

3. The pathwise flux across a boundary. 

Given an open set D we want now to define the flux across dD by 
the path of a diffusion. 

In order to do this we would like to introduce a pathwise analogous 
of the flux as the function 

where Nqj^Ij) (resp. Nqj^^'j)) denotes the number of inward (resp. 
outward) crossing by [0, +oo) 9 t ^-> 7(t) of dD. The problem is that 
the above definition makes no sense since P-a.s. the set {t : Xt & dD} 
has no isolated points and is uncountable. Therefore we are forced to 
proceed in an alternative way. 

Let us observe that if # : 7(t) G dD} < +00 then Ngoi't) is the 
total mass of the random signed measure 

te{s:'y{s)edD} 

where c{t) = +1 if i corresponds to an outward crossing and c{t) — —1 
if t corresponds to an inward crossing. Therefore 

te{s:-f{s)edD} 

where the derivative has to be intended in distributional sense, and 
thus we give the following 

Definition 3.1. Given any open domain D, we define the random 
distribution 
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by 

/• + 00 

(/i^(7), 0) := xoiim 0(0) + / XDilit)) m dt . 



where 6 e Vm = Q 



Supposing now that D is bounded and that Xt is transient, we have 
that P-almost surely the random distribution fin has a compact support 
and so its mass, which we denote by NgD, is well defined. We define 
then the fiux across dD by := K^Nqd). 

By the continuity equation (2.5) such an expectation can be exphcitly 
calculated in terms of the current velocity v (use theorem 7): 

Theorem 3.2. Let {n, P*, X^, P) satijy El and H2, with pt G H\W) 
and Vt G H^{M.^) for any t > 0. For any open bounded domain D , with 
dD a finite union of manifolds, one has 

P+OO P 

^dD = / / p{t, x) v{t, x) ■ n{x) da{x) dt , 

Jo JdD 

where n denotes the outward unit normal vector along dD and a is the 
surface measure. 

4. Diffusion with an asymptotic velocity. 

Since our goal is to define an asymptotic fiux across hypersurfaces, 
we need to impose a condition on the time evolution of the process 

Definition 4.1. We say that the diffusion paths admit an asymptotic 

velocity when 

H3) 

lim - = Woo ^ P-a.s. 

t]oo t 

and moreover jioo <^ A, where A denotes the Lebesgue measure on 
and Poo is the distribution of foo- 

From now on by an open cone Cy, we will mean a set of the form 

{As G : X G S, A > O} , 

where S is an open subset of the unit sphere with (9S a finite union of 
manifolds. 

Remark 4.2. For any open cone Cg, any ball of radius i?, and for 
any diffusion (fi, JF, P^, X^, P) satifying HI, H2 and H3, one has 

li> Xc^nBdXt) = lim XCj^iXt) = Xce(woo) P-a.s. . 

tyoo t]oo 
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See [ini for a two-line proof. Thus hypotheses H3 requires that the 
hmiting velocity is non negligible and such that asymptotically the 
paths have the same direction as their limiting velocity. For example 
the Brownian motion in is transient but it has no limiting velocity 
according to our definition because of the S. L.L.N. Only a Brownian 
motion with drift could satisfy the requirement of our definition. 

A simple but general condition giving the existence of a limiting 
velocity, which is again a finite relative entropy condition (now on the 
full cr-algebra and with respect to P, see the proof below), is given in 
the following (see 

Theorem 4.3. Let (fi, J^, Ft, Xt, P) satify HI and H2. If moreover one 
has: 

r+oo 

E \\b{t,Xt)-Xt/tfdt<+^ to>0, (4.1) 



then: 



lim - Xf = f oo P-a.s. 

t'\oo t 



for some random variable Voo- 

Proof. The condition (4.1) implies, by prop. 2.11, that P -C P 
on a{Xs,tQ < s < +oo), where P is the weak solution of the simple 
stochastic differential equation 

dXt = jXtdt+ dWt . 

with Wt a standard Wiener process and is such that PoX^;^ = PoX^^^ 
Therefore: 

- Xt] =-^Xtdt+- dXt = -Wt . 



and so 



Since 



1 1 /■* 1 ~ 

-Xt = -Xt, + / -dWs 
t to J to s 



Ell - dWs I = / — ds <+oo 



to ^ ) Jto ^'^ 



by Doob's martingale convergence theorem one gets P-a.s. convergence 
of J Xt. The theorem then follows by absolutely continuity. □ 
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Remark 4.4. Under the same hypotheses of Theorem 4.3 it is possible 
to prove (see that the random variable Voo generates the tail a- 
algebra 

T:=f]a{Xs, s>t) . 
t>o 

5. The asymptotic flux across hypersurfaces. 
Let us consider hypersurfaces of the following type: 

where 5*^ is the sphere of radius R. We will define the pathwise flux 
across S/? in the limit when i? | oo. 

Suppose at first that # {t : 7(t) G Sr} < +oo. Since, by H3, t 
7(t) is definitively either in Cs or in C^, if R is sufficiently large one 
has 

te{s : 7(s)eSfl} te{s -. 7(s)eSflU(aCEUB|j)} 

We are therefore lead to give the following 

Definition 5.1. The asymptotic pathwise flux across S is defined by 

:= hm N^, , 

R]oo 

where iV^^ is the total mass of the random distribution — /icEns^- 

The following result shows that the above definition makes sense. 

Theorem 5.2. Let (fi, J^, Pi, Xi, P) satify HI, H2 and H3. Then P- 
almost surely the random distribution Hc^nB'^ has a compact support 
and so its mass Nj^^ is well defined. Moreover one has 

lim iVsfl = XCsi^oc) f-a.s. . 

Rloo 

Proof. Let 

tr{^) :=sup{t>0 : ||7(t)|| < R} . 

By H3, tr < +00 P-a.s. Thus /icnsfj has compact support P-a.s.. 

Let (j)-y G V{R) such that 0^ = 1 on a neighbourhood of [0,rij(7)]. 
By the definition of /icEfiBfj one has 

/• + 00 

(/iCsnB- (7),07) = - XCsnB- (7(0)) - XCs(^oo(7)) / ^^{t) dt 
and the thesis then immediately follows by taking the limit i? t 00. □ 
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The next theorem shows that the definition of asymptotic fiux across 
S by <l>|; := K{N^) is consistent with the result given in theorem 3.2 in 
the case of the fiux across a boundary: 

Theorem 5.3. Let J", Ft, Xt, P) satify HI, H2 and H3 and suppose 
pt e H\R^), vt E H\R^) for anyt>0 and 

e[ \\\/v{t,Xt)fdt <+oo (5.1) 
Jo 

for some T > 0. Then 

i>+oo p 

$^ = E(xc'j,(t;oo)) = hm / / p{t, x) v{t, x) ■ n{x) da{x) dt . 
Proof. Proceeding as is [T^ one has 

p+oo p 

E(xc'e(^oo)) = liiii / / p{t,x) v{t,x) ■ n{x) da{x) dt . 

^Too Jq JT.aUidC'snB'i^) 

The proof is then concluded by proving that: 



-oo 



lim / / p{t,x) v{t,x) ■ n{x) da{x) dt =0. (5.2) 



By the monotone convergence theorem, (5.2) follows from 

/ \p{t,x)v{t,x) ■ n{x)\da{x) dt < +oo (5.3) 

Jdic^nB'^^) 

for some T > 0. Since 

! ! \p{t,xy/^p{t,xY/^v{t,x) ■ n{x)\da{x)dt 

Jo Jo(CEnB^) 

1/2 

< I 11/ p{t,x)da{x) 1 X 

IdiCsnB-j^) 



^ 1/2N 

|2 



pit,x)) \\v{t,x)\\^da{x) I \ dt. 

9(CEns^) 
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by trace estimates on functions in H^(M?) (see e.g.p, chap. 5) one has 

\p{t, x) v{t, x) ■ n{x) I da{x) dt 



Jd{CT.r^B%) 

T ( I- / I- \ 1/2 



< 



/ ( / {p(t,x)dx^ [ \\Vp'/'{t,x)rdx] 
Jo \Jr^ V Jm.^ / 



1/2^ 

x( / p{t,x)\\v{t,x)fdx+ I \\\/{p^^^v){t,x)f dx) ] dt. 



From (2.4) we have u = hence one has Vp^^^ = up'^^'^ and 'V{p^^'^ 



V 



uvp^l'^ + p^l'^SJv. Therefore in order to obtain (5.3) it is sufficient to 
have (5.1) and 

E / (||M(t,Xi)f + \\v{t,X^)f) dt < +00 
Jo 

which is equivalent to (2.2). □ 



6. CONCLUSION 

The present paper introduces the notion of flux across the boundary 
of a bounded open set D in Euclidean space. For regular compact 
domains it is the expectation of the total mass of the (almost surely 
compactly supported) random distribution given by the distributional 
time-derivative of the functional 

where the diffusion Xt, which satisfies the stochastic differential equa- 
tion dXt = bt{Xt) dt + dWt, is transient. If, moreover, the limit 

lim - Xt 

fjoo t 

exists almost surely, then to the notion of asymptotic flux is given a 
sense as well. The flux can be expressed in terms of pt, the density of 
the distribution of Xt, and the current velocity Vt = h — log pt- 
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